A Turing model of pattern formation in atmospheric pressure gas discharges 
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In this letter we propose a Turing model of the formation of patterns of visible light emission 
intensity in atmospheric pressure gas discharges. The electron density and the electron temperature 
take the roles of activator and inhibitor respectively in a two-reactant Turing model, with the 
activator diffusion coefficient being much smaller than that of the inhibitor, and ionization and 
excitation from excited state atoms considered as the dominant reaction processes. The model 
predicts striations in a ID system, which quantitatively agree with experimental results in terms of 
the spatial variation scale lengths. Additionally, the model also predicts changes in the discharge 
structure observed experimentally when input power and gas pressure are varied. 



Plasmas, especially those at atmospheric pressure, of- 
ten exhibit non- uniform and/or non-stationary emis- 
sion intensity of visible light. These phenomena, called 
patterns, are common in non-equilibrium systems [U- 
3 . Patterns including striations, constrictions, zigzags, 
hexagons, clusters, etc. have been observed in gas dis- 
charges since the 1830s [3MlO|. Models based on skin ef- 
fects non- uniform gas heating [HHH] and electron- 
electron collisions JT3] were proposed to explain these 
phenomena, but all these interpretations could not ex- 
plain some of the important characteristics of the pat- 
terns such as the discharge instabilities and the spatial 
periodicity. Simulations based on the kinetic theory suc- 
cessfully reproduced striations [3] and constrictions [13] . 
but the basic mechanism of the pattern formation is not 
revealed due to the many complicated processes included 
in the simulations. 

In 1952, A. M. Turing adopted bifurcation analysis 
of coupled nonlinear reaction-diffusion equations to ex- 
plain pattern formation in reaction-diffusion systems [15] , 
which is now referred to as the Turing model. A typ- 
ical two-reactant Turing model includes an activator- 
inhibitor system, where the diffusion coefficient of the 
inhibitor should be much larger than that of the activa- 
tor [3] . The model was later applied to explain patterns 
in DC discharges [7HT0 ] [16 ] [17] . with current density and 
electric potential chosen as reactants. However, the the- 
ory is limited to the two-layer discharge structure of the 
system, and thus could not explain the patterns reported 
recently in other gas discharge systems, including mi- 
crowave plasmas [18], dielectric barrier discharges [19], 
boundary layer discharges [20] and breakdown processes 
[2"T] . These patterns share similar characteristics, for in- 
stance, when the pressure is decreased or the power is 
raised, the discharge becomes uniform [TM2TJ] . Such con- 
sistency suggests that there should be a common physical 
mechanism for these patterns. 

In this letter we introduce the electron density and the 



electron temperature as the activator and the inhibitor 
respectively into a Turing model to explain pattern for- 
mation in gas discharges. These two reactants are fun- 
damental parameters controlling the various processes in 
plasmas. While this model can be applied to various dis- 
charge systems, we use a particular one at atmospheric 
pressure as an example here to illustrate that the predic- 
tions of the model can indeed be quantitatively consistent 
with experimental results. 

The experimental observations considered here are the 
striations of the visible light emission reported in atmo- 
spheric pressure argon discharges generated with a mi- 
crowave induced plasma source based on a microstrip 
split-ring resonator [18] . where the striations are dis- 
tributed periodically along the gap between the elec- 
trodes. Some experimental images are shown in FIG. 
[2] (a)-(c). This periodic distribution of the visible light 
emission intensity suggests a periodic distribution of the 
excited state atom density, which depends on the electron 
density n and the electron temperature T . So it is implied 
that their distributions should also be non-uniform. To 
describe such non-uniformity, we first consider the elec- 
tron continuity equation: 

dn d 2 n 

Here x denotes the position along the gap, and D a is the 
ambipolar diffusion coefficient, while G n and L n are the 
production and loss rates of the electron density respec- 
tively. At atmospheric pressure, the ionization from the 
excited state atoms dominate the production of the elec- 
tron density, and the collisional-radiative model suggests 
that the production rate is [2"2"] : 

G n = nn 2 e- E ^ T (2) 

Here r\ is a constant and E\ is an effective ioniza- 
tion energy for excited state atoms. In atmospheric 
pressure argon discharges, the electron loss is mainly 
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due to electron recombination with argon dimer ions: 
Ar+ + Ar + Ar -t Ar + Ai 
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[22] . The rates of these reactions can be expressed as 
Ri = r2n{Ar + ) and R 2 = L n respectively, with r 2 be- 
ing a constant. These two processes are also the dom- 
inant production and loss processes of Ar 2 , which are 
locally balanced: R\ = R 2 . In addition, according to 
the collisional-radiative model [22] , the density of Ar + is 
much higher than that of Ar 2 , thus approximately equal 
to the electron density: n(Ar + ) = n. Then we have: 



R 2 = Ri = r 2 n(Ar + ) — r 2 n 



(3) 



The patterns in gas discharges often come with a non- 
uniform distribution of the electron temperature [4], the 
existence of which has been proved possible in the dis- 
charge system discussed here [53] [21] . To describe such a 
no-uniformity, we start with the electron energy conser- 
vation equation: 
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Here p is the pressure, and we assume the energy dis- 
tribution of the electrons to be Maxwellian: p = nT. 
u is the flow velocity and kt = (3/2)nD e is the ther- 
mal conductivity, where D e is the diffusion coefficient 
of the electrons. Ge and Le denote the gain and loss 
rate of the electron energy respectively. In the ambipolar 
regime, D e is 2 ~ 3 orders of magnitude larger than D a . 
So, compared with the heat conduction term, the energy 
flux caused by the flow, that is, terms containing u, can 
be neglected. In the bulk plasma, G E is mainly caused 
by the Ohmic heating |25j . When filaments form, more 
power is coupled to these channels of high conductivity 
(high electron density) than to the neighboring regions 
of lower conductivity (lower eletron density). In a first 
approximation, this non-uniform power deposition can 
be modeled by assuming that the power density is pro- 
portional to the plasma conductivity, and thereby, the 
electron density: 



G E = (3/2)r 3 n 



(5) 



Here r 3 is a constant that depends on the input power. 
As to Le, the excitation of excited state atoms dominates 
the loss of the electron energy [22] : 



L E = (3/2)r 4 n 2 e- B2 / T 



(6) 



Here r^ is another constant, and E 2 is an effective excita- 
tion energy, with E% < E±. Based on the considerations 
above, Q becomes: 
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FIG. 1. Spatio-temporal evolutions of the electron density (a) 
and the electron temperature (b) distributions. The discharge 
condition is 1 atm, 1 W. 



The time derivative of the electron density in the first 
term of the LHS is given by ([I]) and can be neglected be- 
cause D e > D a , r 3 > Tr 2 , and r 4 e~ E2 / T > Tr x e~ El / T 
[22] . We also neglect the third term on the LHS because 
a bifurcation analysis [3[ [TS] shows that this term does 
not affect the periodicity of the striations or the phase 
diagram of discharge transitions (which will be discussed 
later). Neglecting the two terms in ( [7]), dividing the equa- 
tion by n, and substituting ^ and ([3| into equation ([I]), 
we obtain a set of coupled nonlinear reaction-diffusion 
equations for the electron density and the electron tem- 
perature: 



dn 
~di 
dT 
~dt 
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(8) 
(9) 



In equations ^ and ([9]), which describe the gas dis- 
charge system, the electron density promotes its own pro- 
duction via stepwise ionization, while the electron tern- 
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perature self-inhibits: a perturbation that increases the 
electron temperature results in increased energy losses 
due to the increase in excitation processes (third term 
on the RHS of equation Q), driving the electron tem- 
perature down and "inhibiting" the initial perturbation. 
In addition, D e 3> D a . So, the two reactants satisfy the 
conditions of a typical Turing activator-inhibitor system, 
and thus patterns can result from the dynamics of the 
system. 

With the parameters of an atmospheric pressure argon 
discharge |22j . including reaction rates, ionization and 
excitation threshold energies, the equations Q and ([9| 
are solved numerically under uniform initial conditions 
(10 14 cm~ 3 , 1 eV) and fixed boundary conditions (10 13 
cm~ 3 , 0.1 eV) (in fact the boundary and initial conditions 
have little influence on major quantitative characteristics 
of the final striations, including periods and amplitudes). 
The spatio-temporal evolutions of the electron density 
and the electron temperature obtained are shown in FIG. 
[T] The profiles evolve from uniform at the beginning to 
steady periodic structures at the end, which can result in 
striations of visible light emission intensity shown in FIG. 
[2j (B). The opposite phase distributions of the electron 
density and the electron temperature shown in FIG.[T]are 
consistent with simulation results based on the kinetic 
theory [Jj. The magnitude of the spatial period of the 
striations, which has an order of ~ 0.1 mm, agrees with 
the experimental results (FIG. [2] (b)). 

It has been observed in experiments that the number 
of striations increases with the input power (FIG. [2j (b) 
and (c)). This behavior is also predicted by the model: 
with bifurcation analysis [H [3] , the spatial period of the 
striations A can be determined with the parameters in 
(|8) and ((9J: 
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Note that the spatial period of the striations decreases 
when r3, which represents the Ohmic heating term ([5]), 
increases. In other words, the higher the input power, 
the more striations exist. 

In experiments, it has also been observed that when 
the pressure decreases, the discharge structure changes 
from periodic to uniform, i.e. striations disappear (FIG. 
[2j (a)). Such phenomena can also be explained with the 
model: bifurcation analysis [5J [3] also suggests that the 
formation of the striations demands the following two 
conditions: 
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(11) 



(12) 



Here T = (Ei-E 2 )/[log{rir 3 /r2r4)}. When ( 1 1 ) is satis- 
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FIG. 2. Phase diagram: discharge structure transition as 
predicted by the proposed model. Images (a)-(c) show ex- 
perimentally observed distributions of visible light emitted in 
argon microdischarges (scale bar =250 um) sustained between 
coplanar electrodes at 825 MHz at three pressure/power com- 
binations. (A)-(C) are the theoretical results of the emission 
density distributions with the same pressure/power combina- 
tion as (a)-(c) respectively. 



are found to be stable, and thus the striations occur. On 
the other hand, ( [l2| is required for the discharge to be 
stable: when it is not satisfied, either because the input 
power is too small (small r$ in (|5|) or there is a large 
recombination loss rate (large r 2 in ([3])), the discharge is 
difficult to ignite and is not stable. 

Based on the physical mechanisms governing the gain 
and loss terms (5) and ((6), we can assume the 

parameters in the model scale with pressure p and input 
power P as follows: 



p 1 P°,r 2 ~p 2 P°,r 3 ~ p°P 1 ,r i ~ p 1 P°, 



D e ~p- l P°,D a 



P 



LpO 



(13) 



Again, with the parameters of the atmospheric pressure 
argon discharge [22] 



the criteria (11) and (12) can be 
presented in a phase diagram, FIG7]2] In the phase on 
the lower right, where the power is low and the pressure 



fied, periodic distributions instead of the uniform profiles is high, inequality (12) is not satisfied and therefore it is 
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not possible to generate stable plasmas under such con- 
ditions. The phase in the middle is where the discharge 
is spatially periodic, and atmospheric pressure gas dis- 
charge (~ 1 atm, ~ 1 W) is usually in that phase. If 
we lower the pressure (to ~ 0.6 atm) or raise the power 
(to ~ 3 W), we enter the phase on the upper left, where 
inequality (11) does not hold and the discharge is uni- 
form. Despite the simplicity of the proposed Turing 
model, the phase diagram shown in FIG. [2] is in good 
agreement with experimental observations and therefore 
suggests that phase transitions between unstable, peri- 
odic and uniform discharges can be explained in terms 
of the activator and inhibitor roles of the electron energy 
and the electron temperature. 

Even though ^ implies a uniform electric field dis- 
tribution, the existence of periodic structure solutions 
obtained from the model is rather insensitive to this as- 
sumption. In fact, numerical simulation results show that 
the model can still reproduce striations, when ([5| is re- 
placed by Ge oc n 7 , 7 varying from to ~ 1.3. In addi- 
tion, it can be shown that this approximation does not af- 
fect the location of the phase boundary in FIG.[2j across 
which the discharge structure changes from uniform to 
periodic. 

In summary, we have shown that striations and dis- 
charge structure transitions in atmospheric pressure ar- 
gon discharges generated with a microwave source can 
be predicted by the Turing model. The generality of the 
physical principles underpinning this model suggests that 
the proposed model can also be applied to other discharge 
systems, where the patterns share similar characteristics, 
such as, the spatial variation scale lengths decrease as the 
power increases, and the discharge structure also changed 
from spatially periodic to uniform when the power is 
raised or the pressure is decreased [19j[2Q]. In fact, in 
most near-atmospheric pressure gas discharge systems, 
the production and loss of the excited state atoms are 
close to local equilibrium due to their small diffusion co- 
efficients. So based on the collisional-radiative model, 
their densities depend on the local electron density and 
electron temperature. As a result, the following reaction- 
diffusion equations can be used to describe the system: 



dn 
~dt 
dT 

~dt 



= D a V 2 n + f(n,T) 
= D e V 2 T + g(n, T) 



(14) 
(15) 



Where f(n,T) and g(n,T) include the contribution of 
the excited state atoms to the production and loss of 
the electron density and electron energy, respectively. At 
atmospheric pressure, where the electron temperature is 
~ 1 eV, the stepwise ionization/excitation contributes 
mostly to the production of electron density/the loss of 
electron energy, and thus / and g contain nonlinear terms 
with df/dn > and dg/dT < 0. Furthermore, D e > D a 
is also satisfied. So, the electron density and the electron 



temperature act as activator and inhibitor in a classical 
two-reactant Turing model. Therefore similar patterns 
should be expected in systems where the stepwise ion- 
ization/excitation dominate the production of the elec- 
tron density/loss of the electron energy. Extension of the 
model is expected to be able to reproduce the selection 
among different patterns observed in higher dimensional 
systems, or explain more complicated pattern behaviors, 
such as the traveling waves. 
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